Theoretical study of the dynamic structure factor of superfluid "^He 
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We study the dynamic structure factor S{ci, oj) of superfiuid *He at zero temperature in the 
roton momentum region and beyond using field-theoretical Green's function techniques. We 
start from the Gavoret-Nozieres two-particle propagator and introduce the concept of quasi- 
particles. We treat the residual (weak) interaction between quasiparticles as being local in 
coordinate space and weakly energy dependent. Our quasiparticle model explicitly incorporates 
the Bose-Einstein condensate. A complete formula for the dynamic susceptibility, which is re- 
lated to S{q,(jj), is derived. The structure factor is numerically calculated in a self-consistent 
way in the special case of a momentum independent interaction between quasiparticles. Results 
are compared with experiment and other theoretical approaches. 

PACS numbers: 67.40 Db, 61.12.-q 



I. INTRODUCTION 

The dynamic structure factor S{q,uj) observed in in- 
elastic neutron scattering experiments on superfluid *He 
consists of a sharp peak at the Landau phonon-maxon- 
roton energy and a structured continuum at higher ener- 
gies. While there is general agreement that the phonon 
part of the excitation spectrum is due to a collective mode 
of the ^He superfluid (zero sound), the character of the 
observed sharp structure at higher momentum transfers 
has been the subject of considerable debate. This debate 
is reviewed in some detail in the books by GrifflnEJ and 
Glydecl. In particular, the analysis of the temperature de- 
pendence of the dynamic structure factor suggests that 
the roton excitation may be a single particle excitation 
from the condensate, while the excitations at intermedi- 
ate momenta are described as a mixture of colleetive and 
single particle processes (Glyde- Griffin scenarioci). 

As was pointed out by Gavoret and Nozieres (GN 
in the presence of a Bose condensate, the time-ordered 
density-density correlation function ^^{q^uo) of a Bose 
system splits into two parts, 

X'iq.^) = K'-{q.u;)Gi{q,Lo)Kp{q,u;)+xl{q.^)- (1) 

The second term describes density-density fluctua- 
tions of non-condensed particles, while the flrst term ex- 
plicitly contains the single particle Green's functions . 
This term is unique to a Bose-condensed liquid and van- 
ishes in the absence of a condensate. In simple physi- 
cal terms, one may interpret this condensate term as a 
density fluctuation due to the excitation of single atoms 
out of the condensate. According to the Glyde-Griffin 
scenario, this interpretation may be useful in the roton 
region. For small energies and momenta, however, GN 
showed that both terms in Eq. (|^) share the same poles, 
which are due to compressional sound waves. 

Systematic attempts to describe the excitations of su- 
perfluid ^He in the momentum region at and above the 



roton minimum within a phenornenological field-theory 
began with the work of PitaevskiiEl who showed that the 
phonon-maxon-roton curve terminates at twice the roton 
energy due to quasiparticle decay. His elegant ideas were 
developed further by Zawadowski, Ruvalds and Solana 
(ZRS)eI. Starting from a phenomenological Hamiltonian 
written in terms of quasiparticle operators, ZRS found a 
continuum in the one-particle spectral density of states 
lying above the sharp single-quasiparticle peak. In con- 
trast to other approaches the work of ZRS is grounded 
in a field-theoretical analysis similar to the GN theory. 
However, since the density operator in terms of quasipar- 
ticles was not known, ZRS did not calculate the dynamic 
susceptibility. Hastings and HalleyQ attempted to ex- 
tend the ZRS theory to calculate the dynamic structure 
factor by considering only single quasiparticle contribu- 
tions to 5((f, w). Their pessimistic view about ZRS-like 
models may stem from the incomplete density operators 
used in their calculations. A more general exfjcession for 
the density operator was proposed by PistolesiS, however 
without a clear connection to the basic GN theory. As in 
the ZRS approach, the role of the Bose broken symmetry 
remains unclear, because the condensate does not appear 
explicitly within that model. 

The goal of our work is three-fold: (a) to find a quasi- 
particle model for S{q,Lu), which is based on the general 
formalism of GN and applicable in the roton region and 
beyond, (b) to flnd a connection between this model and 
the phenomenological field theory of Pitaevskii and ZRS 
and (c) to calculate S{q, to) consistently (i.e. under con- 
sideration of all terms in Eq. (P) within the model. This 
entails a consistent calculation of the sharp roton peak 
together with the multi-particle continuum. 

In Sec. |l|we discuss the general formalism of GN intro- 
ducing our notation and sign convention. In Sec. [Ill 



we 



introduce the concept of quasiparticles into the micro- 
scopic theory of GN. Treating the interaction between 
quasiparticles as being local in coordinate space, we ar- 
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rive at a model of ZRS-type, but with an explicit consid- 
eration of the Bose condensate and with an expression for 
, which includes both terms in Eq. (^. We study qual- 
itative properties of the model in Sec. Thereafter, 
we present the iteration scheme used in our numerical 
calculations of S{q, cj) (Sec. |v|). The results of the calcu- 
lation are presented in Sec. ^Wtogether with comparisons 
to experiment as well as to other theoretical calculations. 
Finally, in Sec. VI], we discuss our conclusions and give 
an outlook for a further refinement of our calculations. 



II. THEORETICAL FOUNDATION 

In this section we review the GN formalism defining 
our notation |and sign conventions. As is known from 
Belaev's workcl, the presence of a Bose condensate intrct 
duces 'anomalous' propagators. Gavoret and NozieresQ 
systematically extended Beliaev's theory to the two-body 
propagator. For the structure factor, in particular, this 
leads to the first term in Eq. ™). 



A. Response function 

The dynamic structure factor S{q,Lu), which is mea- 
sured in neutron scattering experiments, is given at zero 
temperature by 



S{q,uj) 



1 



(2) 



the ^-function cuts off negative frequency contributions 
and n denotes the density of the Bose fluid. x^{Q:^) is 
the Fourier transform of the time-ordered density-density 
correlation function 

1 r°° 

X^{q,u;)^—J dte-V(9,0, (3) 
which is defined as 

p,p' 

Here, p^{t) = 5^^at__./2(t)ap.+^-/2(i) is the density oper- 
ator in the Heisenberg picture and T the time-ordering 
operator. 

B. Two-particle Green's function 

Following Gavoret-NozieresE, we now consider the two- 
particle Green's function, 



-pi 



The greek indices take the values 1 or 2 and distinguish 
between creation and destruction operators, i.e. a^{t) — 
a^{t) destructs a Boson (^He atom) with momentum q 

at time t and a?.(t) = aLg'(i) is a corresponding creation 
operator. Greek subscripts and superscripts on the left 
hand side of Eq. ^ label incoming and outgoing particles, 
respectively. 

It is convenient to work with the Fourier transform of 
with respect to the time variables: 

(27r)4<S(o,3 +tOi-LO,- lo2)k'^}{p,p', q) = (6) 

Vdi2dt3di4e^^"^*^+"^*^-"=*^-"^*^^i^^^(p,p',g,tO- 

Here,p= (p, e = (^^4 - W3)/2), p' = (p',e' ^ (cj2-t^i)/2) 
are the relative momenta and energies of the particle 
pair in the initial and final states, respectively, and 
q = {q^Lj ^ uj^ + UJ4) is the total momentum and energy 
of that pair. 

It should be obvious from Eqs. (|) and (|) that the 
dynamic susceptibility can be written as 



'Kil{p,p',q), 



(7) 



p,p' 



where we have used the abbreviation J^p = J {2^)'^ I If ■ 
In their diagrammatic analysis of (see Fig. |^), GN 
have shown that the dynamic susceptibility may be sep- 
arated into two terms, a singular and a regular part. 



(8) 



As was already alluded to in the introduction, the sin- 
gular susceptibility Xs (?) arises due to the Bose conden- 
sate and is given by 



xi(g)=A"(g)Gf(g)A^(g). 



(9) 



In Eq. (^) and elsewhere in this paper, summation over 
repeated indices is implied. The singular part of^'^ con- 
tains the Beliaev single particle Green's functionsB G^{q). 
Gonsequently, any structure in G^{q) will show up di- 
rectly in the dynamic susceptibility. This is an interest- 
ing physical feature of a partly condensed Bose liquid: 
excitations of single particles out of the condensate de- 
termine partly the density fluctuations of the liquid. The 
Beliaev Green's functions obey the Dyson-Beliaev equa- 
tion, 

Giiq) = Go^(g) + Go^(g)E«(g)Gf (9) (10) 
with the free Green's functions 



Goliq) 
Goliq) 



Goli-q) = (c. 
Goliq) = 0. 



£0- 



+ p + irj)' 



ill) 



and the self energies E^(g)'s. The free-atom kinetic en- 
ergy is denoted by and the chemical potential by p. 
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FIG. 1: Diagrammatic representation of the two-particle Green's function in energy-momentum space. The solid lines 
represent the single-particle Green's functions G^. The dashed lines indicate condensate particles. F and P stand for the 
four- and three-point vertex function, respectively. The three diagrams without any condensate line are the regular two-body 
propagator F^J given by Eq. (|l^). They are the only contribution to K^J^ in the absence of the condensate. All the other 
diagrams are unique to a Bose-condensed liquid and lead to the singular part of the dynamic susceptibility (see. Eq. (^)). 



The Bose vertex functions A" (q) determine how strongly 
single particle excitations contribute to the dynamic sus- 
ceptibiUty, 

A''{q)^nJ\62c,+6i^)+iY.Poil{p,q)Pte{P,q)- (12) 

p 

Here, uq denotes the density of condensed bosons. These 
functions vanish in the absence of the condensate. The 
interaction vertex (full three point function) P will be 
discussed in more detail below. Here and elsewhere in 
this paper, a boldfaced symbol stands for matrix func- 



tions. The product of two Belaev single particle propa- 
gators G with total momentum q is abbreviated as 

Fol'^{p,q) = G^(-P+ f )G^b+ f ). (13) 

The regular susceptibility is given by 

xS('?) = -*E^2i(P,P','?), (14) 

P:P' 

in terms of the regular two-body propagator 



F^}{p,p', q) - -(2^)^ Foilip, q)S{p p') + Fo^Kp, q)S{p + p') zFo%(p, q)Tf^{p,p', q)F,':^^{p\ q). (15) 



(57 
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The regular two-body propagator is determined by the 
single particle propagator G and the (full) two-body in- 
teraction kernel T. Terms where the condensate does not 
contribute will be called 'regular' in this paper. Xr rep- 
resents the full response function of a Bose liquid in the 
absence of the condensate. 

Following GN we separate out processes with inter- 
mediate two particle states in the interaction kernel F, 
so that it fulfills the following (coupled) Bethe-Salpeter 



equations, 

r^f(p,p',g) = /^f(p,p',g) (i6) 

+ E ^cT(p>p". q)F,fAp", qr'A{p",p', q)- 
p" 

I is the two-particle irreducible interaction kernel and 
contains all diagrams without intermediate two-particle 
states. The full interaction kernel T not only determines 
the structure of the regular susceptibility but also the 
structure of the three-point kernel P and, consequently, 
the structure of the self energy S. Using graphical meth- 
ods one can show, that, in terms of F, the three point 
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function P is given by 

+ i*Ercrb.p'>9)^o;il(y,g)J,"c(y,9), (i7) 



(18) 



and the self energy by 
with 



J and S are again two-particle irreducible functions and 
contain all diagrams without intermediate two-particle 
states. As is obvious from these equations, the vari- 
ous terms G,P, and F, which determine the structure 
of the susceptibility, are not independent, but are linked 
in a complicated way. For instance, any singularity (i.e. 
structure) in T will show up in all terms, and care must 
be taken to treat them consistently. 

In order to put the general formalism outlined above 
to work, the two-particle irreducible functions S, J and 
I must be determined. GN have done this in the hydro- 
dynamic limit q,uj — > 0. But in general, this is a most 
formidable task. Therefore, the formalism is only useful 
if S, J and I can be approximated in a simple way. 



III. QUASIPARTICLE MODEL 

In the following we will introduce the approxima- 
tions appropriate for the momentum regime we are in- 



at all useful, then the residual interaction of the quasipar- 
ticles must be weak. We will exploit this interpretation 
in the following: All the strong atomic interactions are 
contained in S generating the 'bare' quasiparticle spec- 
trum, which, of course, is not at all parabolic as for a free 
He atom. 

Eq. ( |20| ) with the self energy (|l^) precisely corresponds 
to the Dyson equation used by Pitaevskiitl as the basis 
of his evaluation of quasiparticle spectrum near its end- 
point. Moreover, the ZRS phenomenological field the- 
ory is built on a similar basis. Thus one can identify 
Pitaevskii-ZRS quasiparticles with helium atoms renor- 
malized by ^(q). 

In order to further specify the quasiparticle we have to 
define the or alternatively g^. In this work, we will 
take 



IT] 



(22) 



This kind of propagator characterizes a stable quasipar- 
ticle with the 'bare' spectrum cjS and normalization fac- 
tor (residue) A^{q). In general, such a Green's function 
should have a second pole at negative energies, but it 
contributes only a small correction at high energies (if 
kT <C uj'^) and therefore, it can be neglected in our con- 
text. 

Furthermore, we assume that the residues fulfill the 
relation 



Al{q)AUq)-Al{q)Al{q) = 0. 



(23) 



terested in. pWe are guided by the the pioneering work nificantly and leads to a ZRS-like expression for G^iq) 



The same relation holds in the Bogoliubov modeO of a 
dilute weakly interacting Bose gas. It simplifies a formal 
algebraic solution of the Beliaev- Dyson equation (pOt) sig- 



of PitafivskiEl as well as the ZRS phenomenological field 
theoryB. 



A. Quasiparticles 



Aiiq) 



UJ - Lu"^-AJ:*{q) 



with 



Starting from Eq. (y_8|) we transform the Beliaev-Dyson 
matrix equation ( p^ into the form 

Gi{q)^g^^iq)+gliq)^H<l)Gliq) (20) 



AJ:*{q) = A:{q)^*^M- 



(24) 



(25) 



with 



liq) = Goiiq) + Gol{q)mq)gl{q). (21) 



At this point it is convenient to introduce the concept of 
quasiparticles. Eq. (^^ is nothing but an expression for 
Green's functions describing the propagation of particles 
with 'bare' propagators given by the g^{q)'s. Thus, we 
assert that, in fact, g^{q) describe stable quasiparticles, 
which are helium atoms renormalized by the two-particle 
reducible part of the self energy. If this interpretation is 



In order to make contact with the original ZRS expres- 
sions, which contain only one Green's function, we have 
to make some specific choices for the renormalization 
functions A^ in Eq. (|^. Typically, one assumes that the 
ZRS calculation scheme corresponds to keeping only the 
diagonal element G\ of the Beliaev-Green's matrix (see 
e.g. Chap. 10 of Ref. |l|) and takes Al^ A\ = A^^i) and 
A} = 1. This choice is also inspired by the Bogoliubov 
approximation. One may argue that at large wavevec- 
tors the anomalous propagators are much smaller than 
the diagonal ones because the Bose coherence factors be- 
come unimportant as in the Bogoliubov model. However, 
since we do not know the explicit form of these factors 
in our model, we will keep the anomalous Green's func- 
tions in the following. At the beginning of the Sec. Vl, 
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FIG. 2: Two possible 'bare' quasiparticle spectra: the Lan- 
dau spectrum (dashed-dotted line) and the Feynman spec- 
trum (dashed line) compared with the free *He energy /2m 
(solid line). 



As we shall see, this corresponds to an RPA-like approx- 
imation. 

In order to determine the dynamic susceptibility in line 

T ,5 




a 



FIG. 3: Diagrammatic representation of the two-particle 
propag ator W (see Eq. (|^) with all external momenta differ- 
ent from zero. F is the regular two-particle Green's function 
given by Eq. ([l^). 



we will present an alternative choice of the parameters 
A^, which is also consistent with ZRS scheme. 

Finally, we suppose that the 'bare' quasiparticle spec- 
trum is already characterized by a roton minimum and 
the residual interaction does not modify the spectrum 
qualitatively in the momentum region up to the roton 
minimum. This assumption is in the spirit of our ap- 
proach assuming weakly interacting quasiparticles. Since 
an unambiguous determination of the 'bare' spectrum is 
impossible, we will study two modcls-adth different spec- 
tra (see Fig. |^) : (i) Landau spectrumErllJ, which was origi- 
nally chosen in order to explain the specific heat data and 
(ii) Bijl-Feynman spectrumE3, proposed first by Bijl and 
then derived by Feynman. Both spectra are qualitatively 
similar, but differ in the roton and maxon energies. 



with Eq. (M) we first define 



p p,p' 

_ (27) 
From Eq. (p^ ) one finds 

Cil) = fo%iq)+fopJq)+fo$''JqKi{q)fo','liq)- (28) 



Using Eqs. (|2^), (g7|) and (g8|) in conjunction with 
Eq. (17) we can express the reducible part of the self 



(29) 



energy ( [19| ) in terms of the quantity fjjiq), 

^t{q)^W{q)fll{q)J!,{q)- 



B. Response function 

If the concept of quasiparticles is at all useful, the resid- 
ual interaction of the quasiparticles must be weak. Thus, 
the interaction energy in the momentum region of inter- 
est should be much smaller than the kinetic energy of 
quasiparticles and the interaction vertices can be treated 
as being essentially local in coordinate space and weakly 
energy dependent. The Fourier transform of a local inter- 
action vertex is a function of only the total momentum 
transfer, i.e. 



The separation of the dynamic susceptibility (|^) into 
a singular and a regular part is very useful, since it em- 
phasizes the role of the Bose-Einstein condensate, which 
couples the single-quasiparticle propagator G into the 
density-density correlation function x^{q). However, for 
practical calculations it is sometimes more convenient 
to sort the diagrams for K^p{p,p' ,q) differently: Let 

W^'p{p,p' , q) be made up of all diagrams with momenta 
of the external lines different from zero (see Fig. P), 



3{p,q)o,3{q), I{p,p',q)^I{q). 



(26) 



W'^}ip,p\q) = Fil{p,p',q) ^Gi{-p+l)G^p(p+\)Pl^(p,q)G^^{q)P2^{p',q)G'^{-p' + f )G|(p' 



(30) 
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W precisely corresponds to the two-particle. Green's 
functions introduced by Fukushima and IsekiE3. Apart 
from the regular two-particle propagator F it contains a 
term, which rcnormalizes F by taking into account the 
possibility that two quasiparticles interact via the three 
point vertex P and propagate as one quasiparticle in an 
intermediate step. It follows from Eq. ( p^ ) that this sec- 
ond term belongs to Xsil) ™ ®- 

With Eq. ( pO| ) we can write the dynamic susceptibility 
in the form 



X^(g)=no^G^(<j) (31) 

+ [g^ (9) + G?(9)] Jfml\{^) + (9), 
where the function {q) is given by 



(32) 



= /21 (9) + hi\{q)PU'i)G;{<i)P'.H<i)hfi{q)- 

A brief discussion of expression (^ ) is in order: The first 
term describes one-quasiparticle excitations and vanishes 
in the absence of the condensate. The third term cor- 
responds to the direct excitation of two quasiparticles. 
Since in the absence of the condensate reduces to 
the regular two-particle propagator, this term goes over 
into the full response function above T\. The second 



term describes an interference between the one- and two- 
particle channel and it disappears as well in the absence 
of the condensate. Eq. (31) shows how these three terms 
must be combined in order to calculate the susceptibility. 
In earlier literature, the dynamic structure factor was ei- 
ther calculated from the imaginary part of the single par- 
ticle propagator (first term in Eq. (|3l]), see e.g. Ref. |^) 
or from the imaginary part of the two-body propagator 
W21 (e.g. Ref. The interference term was neglected. 
Juge and GriffinO emphasized that all terms may be im- 
portant, however they did not combine them in order to 
calculate S{q, lo). 

It is possible to express the susceptibility entirely in 
terms of g^{q) and f^p{q)- To this end, we transform first 
the expressions ( p8| ) and (|3|) by using of (^ and (|l7|). 



ft}iq) = f^%{q) 



/0/3, 



1 



(33) 
(34) 



Obviously, we have found a set of algebraic equations 
for f^'piq) and w|J((z), which are similar to the Dyson- 
Beliaev equation (j2^) for (g) . We can solve them for- 
mally to obtain 



X^iq) 



a/3 



(l--g(g)J^(g)f(g)J(<r))-i 



9^^{q)+nl^'Y: 



il--{{q)3{q)giq)3^{q)~' 



a/3 



21 



(1 - -g(g)J^(g)f(g)J(g))-i 



9?iq)jfiq)mq) 



(35) 



with fp]{q) given in terms of h%{q) 



fpAq) = 



(1 - -io{q)m- 



a/3 



(/Oa]<(9)+/o7a(9)) 



(57 



per 



(36) 

obvious matrix notation is employed. The denominators 
in (^) clearly show the RPA character of the approxi- 
mation applied here. Note, that all terms in x^iq) have 
the same pole structure. Thus both the sharp peak and 
the continuum observed in S{q,uj) receive contributions 
from all three terms. It may be useful to stress, that 
up to this point, no approximations are involved in the 
derivation of ( ^5| ) apart from Eq. (p6|). However, in order 
to explicitly evaluat e y-^ w e will need the approximations 



for made in Sec. [II A 



We want to end this section by noting that an expres- 
sion similar to Eq. (pW was proposed by PistolesiB as 



an extension of the Zawadowski-Ruvalds-Solana work. 
His Eq. (4) corresponds term by term to our Eq. (|3l]). 
But there are also some differences between both calcu- 
lations. First of all, the Bose-Einstein condensate does 
not appear explicitly in Pistolesi's expression. Instead 
he introduces two fitting parameters. Secondly, he deals 
with only one Green's function and the connection to the 
general formalism of Gavoret-Nozieres remains unclear. 
To evaluate Pistolesi extracts the regular part of the 
two-quasiparticle propagator from experimental data. In 
our calculation (see Sec. |^, we will calculate it within 
an iteration scheme. 
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IV. QUALITATIVE CONSIDERATIONS 



where 



In this section, we will qualitatively study two specific 
kinematical regions, which deserve special attention: the 
endpoint region and quasifree scattering at high energy 
and momentum transfers. 



A. Endpoint of the spectrum 

As has been shown by Pitaevskiii, the renormalized 
single-quasiparticle spectrum of a Bose liquid has an 
endpoint, i.e. at zero temperature undamped excita- 
tions cannot exist at momenta larger than some thresh- 
old value, provided the 'bare' quasiparticle spectrum is 
characterized by a minimum. Pitaevskii has clarified the 
character of the spectrum near its endpoint in a quite 
general way by analyzing the singularities of the single- 
quasiparticle propagator. 

Pitaevskii started from the following expression for the 
self energy (written in our notation): 



(2^) 



d''pP{p,q)G{p)Giq-p)Jip,q) (37) 



which is similar to our Eq. (|19|). The only difference is 
that Pitaevskii did not consider all Beliaev-Green func- 
tions as we do. However, since all these functions have 
the same pole structure (see e.g. Ref. |[), near the sin- 
gularities our expression should lead to the same results 
except for some coefficients and/or regular contributions, 
which are not of interest for the analytical behavior of the 
propagators. Thus Pitaevskii's prediction about the end- 
point of the renormalized spectrum due to the decay of 
quasiparticles holds within our model. 

Pitaevskii distinguished three kinds of decay processes 
into two excitations: (i) phonon creation, (ii) decay into 
excitations with finite momenta propagating in the same 
direction with the same velocity and (iii) decay into two 
rotons. 

Since we are only interested in the momentum region 
10 < <f < 40 nm^^ and do not take phonons into account, 
the emission of two rotons is the only possible decay chan- 
nel. Thus, the spectrum should have an endpoint of the 
third kind, which is characterized by a logarithmic sin- 
gularity. This endpoint leads to a threshold at twice the 
roton energy in the imaginary part of the self energy, i.e. 
to 



Iui^*'i{q,Lu < 2 A) = 0. 



(38) 



Here, A denotes the roton (minimum) energy. 

Using Eq. ( |38| ) in conjunction with Eq. ( p4| ) we find 
that the imaginary part of the one-quasiparticle Green's 
function consists of a sharp peak below and a continuum 
above twice the roton energy. 



ImG^(g) = -nA^{q)6iLO-Lu°ff-AReE*{q))0{2A-i 



-ImM^iq), 



ImA/f (g) = 



< 2A 



(40) 

To proceed further we use S{f{x)) = J2i ^{^ ~ Xi)/ f'{xi) 
with f{xi) = and assume that the equation 

UJ-uj"^~ARe^*{q) = (41) 

has only one solution ui^ below 2A. We then arrive at 

ImGf (g) — — TrZf (g)J(cj-w ^-)6l(2A-a;)-HlmM£((j) (42) 

with a renormalization factor given by 



Ai{q) 



II 



(43) 



We see that in fact, the sharp component of the one- 
particle propagator can only exist below the threshold 
energy 2A as predicted by Pitaevskii. Furthermore, the 
only effect of quasiparticle renormalization below the 
threshold is a modification of the excitation energy and 
peak strength due to the real part of the self energy. 
There is no mechanism in a model without phonons, 
which could change the width of the peak below the 
threshold. 

It is interesting to analyze the behaviour of the peak 
strength near the endpoint. To this end we first need 
to determine some properties of the real part of the self 
energy AReTi*{q) in the vicinity of the threshold energy. 
Since the real and imaginary part of the self energy are 
related by 



^ReE*(g,w) 



1 



d.^1^2^1^, (44) 



one can show that the threshold in AImI]*(q) leads 
to a logarithmic singularity in ^ReE*(g), schematically 
shown in Fig. IJ. At the threshold energy the function is 
finite but its left and right derivatives are infinite; 



aylReE*(g,u;) 



if 



2A. 



(45) 



It follows immediately from Eq. ([45|), that the peak 
strength vanishes at the endpoint energy. 



(46) 



Recall that the dynamic susceptibility x^iq) shares poles 
with the single-particle Green's function. Thus we can 
expect, that the strength of the sharp peak in the imag- 
inary part of x'" {q) (i.e. in the dynamic structure S{q)) 
vanishes at the threshold as well. Such a behavior is in 
qualitative agreement with the standard interpretation 
of the experimental data. 



(39) 
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FIG. 4: Logarithmic singularity in AReT.*{q) due to the 
threshold in AIm'S*{q) at a; = 2A. At the threshold energy, 
the function is finite, but its left and right derivatives go to 
infinity. 



B. High- momentum scattering 

In the high-momentum region, S{q, u) can be increas- 
ingly well described within the impulse approximation. 
In the following we will discuss the mechanism leading to 
this approximation within the framework of our model. 
For a more general consideration we refer to Refs. |l|, ^ 
and 

We start by noting that if the momentum is high 



enough the kinetic energy of the quasiparticles must be 
very large relative to the potential energy of their interac- 
tions. Thus interaction effects (i.e. all terms containing 
I and J) can be neglected and we obtain from Eq. (|35| ) a 
simplified expression for the dynamic susceptibility, 



a0 



af3 



UJ 



(47) 



with 



/oaf te^) 



(2^)4 



dpy deg„(p,e)5^((f-p,w - e) 

_ (48) 
and given by Eq. (P3). Here, we have only two con- 
tributions to ■ The first term proportional to uq in- 
troduces single-quasiparticle excitations into the density 
fluctuations and produces a sharp peak in S{q,uj). The 
regular part Xr leads to a continuum in S{q, ui). It is due 
to simultaneous excitation of two non-interacting parti- 
cles. Unlike in Eq. (p5|), the terms contributing to 
do not have the same pole structure. In Ref. |l|, a result 
similar to Eq. (4^, however with different 'bare' Green's 
functions, has been called an improved form of the Bo- 
goliubov approximation. 

In Eq. (Eq), we carry out the integration over e ana- 
lytically using the theorem of residues. After that the 
dynamic susceptibility can be written as 



X (g,c.) = no>^<7f(9,-)+y ^ ^_^o__^o + ,^ ■ (49) 



Finally, we find that 



5(g,a;) = ^5]A^(g)<5(c.-c.a) + i / -^[Ajmliq - p) + AjmUq - p)]Siu^ - 4 - colp). (50) 

If the dominant contribution to the integration over p in Eq. ( |50| ) is from momenta much less than q, we can use 
J d^p ~ d'^p, where D = {p : \p\ I'll}- this region, we have + — ^-p- large enough momenta, 
renormalization effects coming from the irreducible part of the self energy S((f, cj) are negligible in comparison with 
the kinetic energy of the particles and we can use the approximation 

^ 2m 2m m 

Here, m is the free ''He mass. We arrive at 

5™^(<f,c.)^^^A^(g)5(a--|^) + i / ^[Alip)Aliq-p)+Alip)Aliq-p)]Sic.-^ + P^). (52) 
n ^-^ 2m n J (27r)-' 2m m 



It follows from the second line of Eq. (p3), that in the high momentum region the continmmr part of S{q,LLj) 
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FIG. 5: Experimental results for ^((f, (.<j)Ei at g = 32 nm~^. 
The solid vertical line shows the position of the peak inter- 
preted as a quasifree peak. The dashed line indicates the 
corresponding free-atom energy. 



consists of a Doppler-broadened peak centered at the en- 
ergy q'^/2m. In other words, at very high momenta the 
continuum is dominated by the free atom properties of 
the system. 

Eq. (52) is in line with the Green's function reformula- 
tion of the impulse approximation presented by GrifhnEJ. 
Some differences are due to specific assumptions about 
the form of made within our approach. 

For wavevectors |g| > 32 nm"-'^, one observes_a peak 



centered close to the free-atom energy (see Fig. |). The 
larger the momentum transfer the smaller the shift be- 
tween the peak position and the free-atom energy. There- 
fore, in ^He , the impulse approximation becomes valid 
already at moderately large momentum transfers. Since 
a small shift in the energy of the observed 'quasifree' 
peak indicates that the difference between uj'i and /2m 
is small as well, the approximation provides a criterion 
for choosing a spectrum of 'bare' quasiparticles. We will 
exploit this in Sec. VI B| . 



V. 



PREPARATION OF THE NUMERICAL 
ANALYSIS 



The formulas of Sec. Ill provide a set of equations 



enabling a numerical calculation of the dynamic suscep- 
tibility. In this section we will briefly discuss some details 
of its solution. 



G 



FIG. 6: Iteration scheme. See text for more details. 



A. Iteration scheme 

We begin by writing down the complete set of equa- 
tions for the calculation of in a well-ordered form, 



Q/3 



1/2 



Y.Gi{q)jf{q)ff^{q)+wl\{q), 





(1- 


^lliq) = 


(1- 


Ciq) - 


(1- 


foa}{q) = - 


i 

f27r)4 



9^p{q), 



rT\-l 



(foil + foil), (56) 



(53) 

1 

(54) 
(55) 



Q/3 



We see from the above equations that the functions 
foi^piq) are all we need in order to calculate the dynamic 
susceptibility x^{q)- However, determination of these 
functions is not a trivial task, since they are contained 
implicitly on the right hand side of Eq. (|57|). We will 
solve these integral equations self-consistently via the fol- 
lowing iteration procedure: We start with = (i.e. 



*iiq) = 0), and calculate fo%. 





The result is used to 



evaluate /^^ from which new G^ are obtained. With 

these we again calculate /o^^ and the calculation is 
repeated until self-consistency is established. The proce- 
dure is illustrated schematically in the diagram in Fig. |^. 
Note that technically the iteration cycle has-similarities 
with a calculation made by Gotze and LiickalJ who con- 
sidered density excitations in superfluid ^He using Mori's 
formalism. 

Obviously, we need the 'bare' quasiparticle propaga- 
tors g^iq) (i-e. (jj'i and the strengths A^{q)) and the 
interactions l{q) and 3{q) in order to calculate ■ For 
the sake of simplicity, we will treat them as independent 
parameters of the model. Neglecting the connection be- 
tween the input quantities has some consequences. For 
instance, if one changes only the 'bare' spectrum uj'i and 
keeps the interactions fixed, one arrives effectively at a 
different physical system. In other words, the result of 
the iteration depends on the 'bare' spectrum. That is the 
reason why we need some criteria like high-momentum 
scattering (see Sec. [VB) in order to choose a suitable 



input spectrum for our calculation. 
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We will use the resulting one-quasiparticle spectrum 
Wq-, which is given by Eq. ([4l|), as a consistency check: a 
solution is consistent, if after some number of iteration 
steps the renormalized spectrum does not change any 
more. 



tically. 

According to Eqs. (| 
into three terms: 



fo%{q) = ft{q)+fri<l) + fr''i'l)- 



and (HI) we can split fo%{q) 

1.1, \ , pl.rni 



(58) 



B. Calculation of /o^XC'z) 



In order to facilitate our iteration scheme we want to 
carry out the integration over energy in Eq. (^7|) analy- 



For the sake of simplicity we have suppressed the greek 
indices on the right hand side of the above expression. 
We find (recall that p — {p, e) and q — {q, uj)) 



J 



/o''(9) 



d^p f de Ai{p)e{2A-e) 



{2n) 



AVq-p)e{2A-LU + e) 



2ti e-uj^p- AKeT,*{p) + irj lu - e - 



■col - 

q-p 



AReJ:*{q-p) + IT]' 



(59) 



r 



for the single-quasiparticle-single-quasiparticle term, 



/o'™(9) 



_fp_ 

(2^ 



de 
2^ 



Ai{p)A'giq-p)ei2A~e) 
e-uj'l- AReJ:*{p) 



irj 



Mi{p)A}{q^p)Q{2A~u;^ 
-e-Lol_p- ARc^*{q-p) 



(60) 



for the singlc-quasiparticle-continuum contribution and 
finally 



Jq (9) = « 



d^p 

(2^ 



de 
2^ 



M^ip)M;iq-p) (61) 



for the continuum-continuum part coming from the prod- 
uct of two singlc-quasiparticle propagators. Here, we 
have introduced the auxiliary function 



Mf(g)-G^(<z)e(^-2A). 



(62) 



Note, that the imaginary part of Mj!^ is given by Eq. (4C ). 

In order to carry out the energy integrals in Eqs. (59)- 
( |6l| ) we first make an analytic continuation of the inte- 
grands into the region of imaginary energies and then 
perform a contour integration. We are dealing with re- 
tarded propagators, so the contour should run along the 
real energy axis and be closed by a large arc in the lower 
half plane. Since the Green's functions behaves like u!~^ 
for \lu\ 00, Jordan's lemma ensures that the contribu- 
tion from the arc at infinity vanishes and we can evaluate 
the integral along the real axis by using of the theorem 
of residues. 

Apart from /"''"(?), which is analytical in the whole 
energy plane, all other terms in /o^^ have single poles in 
the lower half plane. After calculating the corresponding 
residues we get for the imaginary part of /o^^ 



Woa]j(9) 



d^p 



'vr J j^{Z'^{p)Zj{q-^ p)5{^ ~ ~ ^,--p)e(2A - ujp)Q{2A - {u ~ up)) 

u^^~)Zl[q- p)e(2A - uj^^p)} (63) 
I 



p)ImMj(^- p, - u;p-)e(2A - Wp-) -I- ImAf^(p, < 



with and up defined in Sec 



IV A 



Obviously, the 

continuum-continuum contribution to Im/o^^j vanishes 
after carrying out the energy integrals. 

We could also calculate the real part of /o^^ in an 
analogous way. However, since the function is analytic off 
the real axis, its real and imaginary parts are connected 



through the causality relation 



de 



(64) 



where P denotes a Cauc hy p rincipal value. Since a calcu- 
lation of the expression ( |64| ) requires less numerical effort 
than the direct method based on Eqs. (^9|)-(|6ll), we will 
use it in our numerics. 
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VI. NUMERICAL RESULTS 



Our calculation scheme contains several parameters: 
the condensate density Uq, the matrix functions I, P, 
and A, and the free quasiparticle spectrum uj^. For the 
condensate fraction we take no = O.ln, where n is the 
density of the liquid. This assumption is in agreement 
with the-condensate density extracted from experimen- 
tal datata a nd s upported by Path Integral Monte Carlo 
calculationalBEJ. |_. 

In a recent papered, we have presented first numerical 
results for S{q, uj) obtained in a simplified version of the 
calculation with a momentum independent self energy 
of the quasiparticles and constant interactions between 
them. Since such an approximation led to some inconsis- 
tencies in the calculated spectra, in the present work we 
have removed this restriction for the self energy. How- 
ever, for the sake of simplicity, we still treat the inter- 
action functions as momentum and index independent, 
i.e. 



(65) 



Here, (74 and 33 are phenomenological constants. Sim- 
ilarly, we assume that the renormalization functions 
A^((7)'s are momentum independent and equal in all 
channels (i.e. A^{q) — A for a, /3 = 1, 2). Obviously, this 
is not in line with the standard choices made to make a 
contact with ZRS-like mod els fr om the microscopic the- 
ory (see discussion in Sec. HI A). However, under these 
assumptions, the elements of the Beliaev-Green's matrix 
function are identical. If we now carry out summations 
over greek indices in our expressions, we are able to ob- 
tain a scheme of ZRS type. 

In each iteration step, we have evaluated the dynamic 
structure factor fo r 56 |(7|-values and 501 w-values. As 
mentioned in Sec. V A , we have used the resulting spec- 
trum Lo^ as a consistency check: a solution was consistent, 
if after a certain number of iteration steps the renormal- 
ized spectrum did not change any more. Depending on 
the input parameters 10-50 steps were necessary to ob- 
tain a consistent solution. 



A. Pitaevskii singularities 

Let us first consider some general features of the results 
within a specific example with the Landau spectrum as 
the 'bare' quasiparticle spectrum wS. 

Results for the renormalized single-quasiparticle spec- 
trum Wq-, which is defined by Eq. ([4l|), are shown in 
Fig. (^. is characterized by an maximum at the en- 
ergy Sa/ — 9.6K, a minimum at A;? — 7. OK and an 
endpoint at (2(f7j, 2Afl), where is the wavevector of 
the roton minimum. At higher momentum transfers the 
renormalized spectrum differ qualitatively from the 'bare' 
spectrum. It bends toward 21^r due to quasiparticle de- 
cay. Since the density of states is high near the local 



extrema of the spectrum, we expect to observe a struc- 
ture in the continuum part of S{q, oj) , whose origin is as- 
sociated with the roton-roton, roton-maxon and maxon- 
maxon pair excitations. They are as well due to quasi- 
particle decay as discussed by Pitaevskiitl. Naively, one 
could expect to also observe features stemming from the 
plateau near the endpoint. But excitations associated to 
this part of the dispersion curve have a vanishing weight 
(the dashed-dotted line in Fig. ^ and corresponding 
singularities will not appear in the pair excitation spec- 
trum. 

In Fig. |[ we show results of the same calculation for 
S{q,uj) at q — 23 nm~^. The dynamic structure factor 
consists of a sharp peak lying at the renormalized single- 
quasiparticle energy lo^ and a continuum, which is not 
featureless at all. As expected, we observe peaks at en- 
ergies 2 Aft., Aji + Am and 2 Am, respectively. Moreover, 
there is structure, which we cannot assign to particular 
regions of the single-quasiparticle excitation spectrum. 
We believe this to be a numerical defect. However, it 
has not been eliminated, since it is not important for the 
following discussion. 

As mentioned above, we do not take phonons into ac- 
count in our calculation. This results in a gap between 
the sharp peak and the continuum in S{q, ui) . For the 
same reason, there is no mechanism in the model, which 
would change the width of the sharp peak below the de- 
cay threshold. 



B. Landau vs. Feynman spectrum 

In this section, we want to compare results for S{q,Lu), 
which differ only in the 'bare' spectrum uj'i from each 
other (Fig. (||) and (|l0|)). All other input parameters are 
the same. They were chosen to reproduce the observed 
dynamic structure factor at q = 23 nm^^ in the case 
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FIG. 7: The renormalized excitation spectrum as calculated 
from Eq. (^ij) (solid curve) and the resulting strength Z{q) of 
the sharp peak in ImG(g) (dashed-dotted line). 
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Fig. ^ suggests that the Landau spectrum is a better 
choice for our calculation. However, as we know from 



'2 



FIG. 8: S{<],uj) at g = 23 nm~^ calculated using the Landau 
free quasiparticle spectrum. The structure observed in the 
continuum is partly due to the roton-roton, roton-maxon and 
maxon-maxon pair excitations. They correspond to the van 
Hove singularities in the density of states function. 



of the Landau 'bare' quasiparticle spectrum as well as 
possible. The values are as follows: 



53 = 0.5 Knm^/^ 



9i 



-0.04 Knm^ ^ = 0.25. (66) 



To include experimental resolution we have convoluted 
the calculated S{q,uj) with a Gaussian function. The 
standard deviation of the Gaussian was taken to be 
2.46 K in agreement withJiie experimental resolution 
obtained by Fak and Bossytj. 

In Fig. ^, results for S{q,uj) at g = 23 nm~^ are shown 
in comparison to experiment. We see that, indeed, in the 
case of the Landau 'bare' spectrum (left figure), the cal- 
culated results describe the two peak structure present 
in experimental data well. The lower energy peak cor- 
responds to single-quasiparticle excitations. The higher 
energy structure comes from the Pitaevskii singularities 
discussed in the previous section. They have merged into 
a broad peak due to the convolution with a Gaussian. 
The theoretical and experimental strengths of the peaks 
agree to a reasonable accuracy. However, the calculated 
peaks appear at somewhat different energies as compared 
to the experimental data. 

The results in case of the Feynman spectrum are much 
worse. We still observe a two peak structure, but the 
peaks are shifted to higher energies. Moreover, the 
strength of the single-quasiparticle peak is too small in 
relation to that of the second one. The shift in the ener- 
gies indicates that the input parameters were too small 
to force the spectrum down to the experimental one. Ac- 
tually, one can try to find some parameters, which lead 
to better results with the Feynman spectrum. But the 
stronger the vertex functions, the smaller the difference 
between the potential and kinetic energy of quasiparti- 
cles and the worse the approximation about locality of 
the vertices (p6|). 



Sec. IV B, it cannot describe the quasifree peak, which 
becomes dominant for g > 32 nm~^, very well. The 
results at g = 32 nm~^ are presented in Fig. |l^. In 
case of the Landau spectrum (left figure), we observe two 
peaks in the calculated S{q,u!). The first one is due to 
Pitaevskii's singularites. Its position agrees with the ex- 
perimental one to a good accuracy, however its strength 
does not fit the measured one any more. We will come 
back to this feature later. The continuum is dominated 
by a peak centered at a; « 35 K. It turns out, that this is 
the Landau energy of an excitation corresponding to the 
wavevector q = 32 nm~^ and the peak can be interpreted 
as quasifree peak. However, it appears at energy 25 K 
smaller than the observed one. 

In case of the Feynman spectrum we still observe the 
shift of the whole structure to higher energies. But as 
expected, now the quasifree peak is qualitatively better 
described. Note, that now the shift is not the same for 
both peaks. The position of the quasifree peak is deter- 
mined by the 'bare' spectrum, while the first peak (due 
to Pitaevskii singularities) comes from the renormalized 
one. So, its position should depend on other input pa- 
rameters, i.e. vertex and renormalization functions. 

As we see, neither the Landau spectrum nor the Feyn- 
man spectrum appears to be a good candidate for the 
'bare' quasiparticle dispersion relation, if one wants to 
describe S{q,u!) in a wide region of momentum transfers. 
While the Landau spectrum is better in the kinematic 
region, where the Pitaevskii's singularities play the dom- 
inant role, the Feynman spectrum is more suited for the 
quasifree region. Since our goal is a description of both 
regions, we decided to modify the 'bare' spectrum. In the 
following we use a phenomenological spectrum given by 
the dashed-dotted curve in Fig. |ll|. The phonon-maxon- 
roton part of the new co'l agrees with the Landau spec- 
trum (solid line), but it approaches the Feynman energy 
(dashed line) much faster above the roton minimum. 

Apart from the Pitaevskii's singularities and the quasi- 
free peak, in case of the Feynman spectrum there is an 
indication for a third peak lying between them (Fig. ^0|). 
It corresponds to the third experimental peak observed 
at w w 39K. But its strength is much to small compared 
with experiment. We will discuss this issue in more detail 
below. 



C. S{q,uj) with 'new' quasiparticles 

The effects of the modified 'bare' quasiparticle spec- 
trum on S{q, iv) are shown in Fig. ^ In this calculation, 
we used 

33 = 0.3 Knm^/^ = 0.19 Knm^ A = 0.4 (67) 

and the Gaussian resolution function discussed in the pre- 
vious section. As expected, the results at g = 32 nm^^ 
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FIG. 9: S{q,uj) at q — 23 nm ^ calculated with different 'bare' spectra: the Landau spectrum (left figure) and the Feynman 
one (right figure). Interactions and g4 and renormalization function A are the same in both cases: — 0.5 Knm^/^ 
P4 = -0.04 Knm^ and A = 0.25. 
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FIG. 10: S{q,ijj) at g = 32 nm ^. See caption of Fig. U for more details. 



are now much better than in the case of the Landau spec- 
trum. First of all, we see that the quasifree peak is sitting 
at the right position. Moreover, weights of the dominant 
peaks agree to a better accuracy with experiment. There 
is also evidence for the third peak observed in experi- 
ment. But its intensity is still much to small. However, 
the results at g = 23 nm^^ and q = 28 nm^^ are not 
satisfactory. In both cases the qualitative structure of 
S{q,uj) is still reasonable, but the agreement with the 
experiment is now much worse than at q = 32 nm~^. 
From the analysis of the results we are led to the con- 
clusion, that it is difficult to describe S{q, lo) in a wide 
momentum region in the framework of ZRS-like models 
with momentum independent input parameters. We will 
remove this restriction consistently during further refine- 
ment of the calculation. 

If the above conclusion was true, we should be able 
to find another set of parameters, which leads to better 
results for S{q,uj) at g = 23 nm^^ or q ~ 28 nm"-'^. 



Indeed, we have found such parameters, which lead to 
reasonable results in the region of smaller momenta. In 
Fig. n3, we present S{q,Lj) calculated with 



33 — 0.5 Knm 



3/2 



54 



-0.08 Knm'' and A = 0.265. 

(68) 

Now, S{q,uj) at q — 23 nm ^ is described quite well. 
However, if we compare Fig. |lj with the results of the 
previous section obtained with the Landau spectrum, 
we see, that the modification of the spectrum results in 
an unreasonable plateau at higher energies. This again 
points to a possible momentum dependence of the inter- 
action parameters. Although somewhat worse than at 
q = 23 nm~^, the calculated spectrum at g = 28 nm^^ 
is still reasonable compared to the experiment in this 
special case of input parameters. And the results at 
q = 32 nm~^, which are now worse than in the previous 
case, seem to confirm our conclusion that input parame- 
ters must be momentum dependent. 

From the input parameters reported in ( |67| ) and (|6S 
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FIG. 11: Our 'free' quasiparticle spectrum (dotted-dashed 
Une): we modify the Landau curve (solid line) in the region 
above the roton minimum such that it approaches Feynman 
spectrum (dashed line) much faster. 



an interesting feature of the irreducible four-point vertex 
function 34 follows. We have achieved the best fit to the 
data at intermediate momenta for a negative 54. But at 
high momenta, 34 has to be positive in order to describe 
S{q,u}) properly. Thus, the effective interaction between 
quasiparticles changes its character from an attraction at 
intermediate momenta to a repulsion at higher ones-. We 
will come back to this feature in forthcoming workEj. 

As we know from Eq. (^l[), there are three contri- 
butions to S{q,u!), which come from the single-quasi- 
particle, interference and two-quasiparticle terms. We 
will refer to them as S^^\ g(mt) ^^^^ 5'^^-', respectively. 
The energy integrated strengths of these contributions 
are reported in Table |. There are two important features 
shown in the table. First, reduces the strength of 

the sharp single-quasiparticle peak at g = 23 nm^^ and 
of the continuum at g = 32 nm~^. Secondly, and more 
important, in all cases the structure factor is dominated 
by the two-particle contribution S^^\ Thus the expres- 



S{q,uj) oc ImG'((f, w). 



(69) 



which was used in the literature (e.g. Ref. ^ in order 
to calculate S{q, uj) within ZSR-type models, cannot de- 
scribe the data in the momentum region under consider- 
ation, since the most relevant term is neglected. 



D. Comparisons with other calculations 

In this section, we want to compare our results with 
those obtained by other methods. We will limit ourselves 
to the work of Q'jp^ and LuckecZI and of Manousakis and 
Phandaripandeoo. All theoretical results are compared 
to the experimental data of Fak and Bossyt£l. 

Gotze and Liicke presented a detailed analysis of 
S{q,uj) for supcrfluid ^Hc at T = within the memory 
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FIG. 12: S{q, ui) calculated from the modified Landau spec- 
trum as 'bare' quasiparticle spectrum. The input parameters 
were chosen to give a reasonable fit to the experimental data 
at (7 = 32 nm~^ gs = 0.3 Knm^'''^, 54 = 0.19 Knm^ and 
A = 0.4. 



function formalism. Their results at two selected mo- 
menta are presented in Fig. |l^. In order to make compar- 
isons with experiment we have convoluted the calculated 
spectr a with the Gaussian resolution function discussed 
in Sec |VIB|. 

In Fig. |15|, we show S{q,uj) calculated by Manousakis 
and Pandharipande, who used correlated-basis- function 
methods at T = 0. Again, the theoretical spectra are 
convoluted with the above mentioned Gaussian resolu- 
tion function. 

In order to facilitate comparison of the different meth- 
ods our best results for the dynamic structure factor are 
plotted again in Fig. 



16, Recall, that each plot corre- 



sponds to different parameter sets in order to simulate 
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FIG. 13: S{q, ui) calculated from the modified Landau spec- 
trum, but with different interaction parameters: gs = 
0.5 Knm^/^, = -0.08 Knm^ and A = 0.265. They were 
chosen to give a reasonable fit to the experimental data at 
g = 23 nm~^. 

an effective momentum dependence of our parameters. 

It follows from the Figs, p^^l^ that all three methods 
are significantly missing strength in the energy region 
between 30 - 50 K at momenta around 30 nm~^. This 
fact may indicate that there is a lack of understanding of 
the underlying physics. 

VII. CONCLUSIONS 

In this paper, we presented a calculation of the dy- 
namic structure factor of superfluid ^He in the interme- 
diate and high momentum transfer region within a model 
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^ ( i Ti t ) 


-0.23 


0.23 


0.01 


-0.11 




0.77 


1.03 


0.87 


0.57 



TABLE I: Energy integrated contributions to S{q, u) from the 
single-quasiparticle, interference and two-quasiparticle terms 
at different momentum transfers. The contributions at g = 
23 nm~^ and g = 28 nm~^ have been obtained with gs = 
0.5 Knm^/2, 54 = -0.08 Knm^ and A = 0.265. For q = 
32 nm"\ we used 53 = 0.3 Knm^/^ 54 = 0.19 Knm^ and 
A — 0.4. The dynamic structure factor at g = 28 nm~^ and 
g = 32 nm~^ contains only a structured continuum. At these 
momenta, there is no sharp peak in ^(g, cj) , since the strength 
given by Eq. (^) vanishes. 

based on the Gavoret-NoziereJi microscopic theory. Af- 
ter the introduction of quasiparticles, we obtained an 
RPA-like expression for the density-density correlation 
function in a model, which has similarities to the phe- 
nomenalogical field theory of Zawadowski, Ruvalds and 
SolanaB, but treats the condensate explicitly. We evalu- 
ated S{q,uj) numerically for the special case of momen- 
tum independent interactions between the quasiparticles. 
Our results suggest the following conclusions: 

(i) The quasiparticles, that were introduced here, can 
be interpreted as helium atoms renormalized by 
the two-particle irreducible part of the self energy. 
Quasifree scattering may be used to determine the 
high energy part of the 'bare' quasiparticle spec- 
trum appropriately. 

(ii) S{q,u!) consists of three terms: one- and two- 
quasiparticle excitations and a interference term, 
S^-^\ 5'^' and S'(*"*\ respectively. All terms have 
the same pole structure and appear to be equally 
important in the momentum region at and above 
the roton minimum. Thus the expression S{q,uj) cx 
ImG((f, w) often used in literature in order to cal- 
culate the dynamic structure factor in ZRS-type 
models neglects important terms. In the absence 
of the condensate only the S'^^-' term remains. 

(iii) A model which employs momentum independent 
interactions cannot account quantitatively for the 
neutron scattering data in a wide momentum re- 
gion. 

(iv) The two-quasiparticle interaction 174 should be at- 
tractive at intermediate momenta and repulsive at 
higher momentum transfers. 

(v) As is seen from the comparison with other calcula- 
tions, our model provides an alternative description 
of the experimental data. However, no method de- 
scribes qualitatively the data at |(f| > 30 nm^^. 
This may indicate a lack of understanding of the 
underlying physics. 
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FIG. 14: Calculation of S{q, u) by Gotze and LiickeQ using the memory function formalism at two selected momenta c ompa red 
to experiment. The calculated spectra have been convoluted with a Gaussian resolution function as discussed in Sec. 
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FIG. 15: Calculation of S{q,uj) by Manousakis and PhandaripandeE3Efl using correlated basis function methods. See caption 
of Fig. |l^ for more details. 
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FIG. 16: Our best results for S{q,uj) at two selected momenta. The case g = 23 nm ^ corresponds to gs — 0.5 Knm"^/^, 
54 = -0.08 Knm^ and A = 0.265. For q = 32 nm"\ we have used ga = 0.3 Knm^/^, g^ = 0.19 Knm^ and A = 0.4. 



We now turn to a somewhat more detailed discussion 
of the point At present, we do not know very well 33 
and (?4 as functions of q. Thus for the sake of simplicity 
we assumed, that they are weakly momentum dependent 
and can be approximated by constants. The same holds 
for the residues A^. For the chosen parameters, we eval- 
uated S{q,LLj) at 56 values of \q\. Our results show that 



it is difficult to describe S{q, lo) in a wide momentum 
transfer region within a model with constant parameters 
53, 54 and A. To avoid this problem, we have 'simu- 
lated' a momentum dependence by calculating the struc- 
ture factor with different parameter sets. For each set, 
we picked out from the 56 spectra those with reasonably 
good agreement with experiment. In this way, we got 
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sets of parameters, that led to good results in different 
momentum intervals. 

Clearly, we have made progress towards the develop- 
ment of a quasiparticle model, which treats condensate 
and non-condensate terms in S{q, to) on the same foot- 
ing. However, the numerical calculation needs further 
refinement which will be addressed in future work. 
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